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ABSTRACT 

We revisit theoretical and observational constraints on geometrically-thin disk 
accretion in Sagittarius A* (Sgr A* ). We show that the combined effects of mass 
outflows and electron energization in the hot part of the accretion flow can deflate 
the inflowing gas from a geometrically-thick structure. This allows the gas to cool 
and even thermalize on an inflow timescale. As a result, a compact, relatively 
cool disk may form at small radii. We show that magnetic coupling between 
the relativistic disk and a steady-state jet results in a disk that is less luminous 
than a standard relativistic disk accreting at the same rate. This relaxes the 
observational constraints on thin-disk accretion in Sgr A* (and by implication, 
other Low- Luminosity Active Galactic Nulcei, LLAGN). We find typical cold gas 
accretion rates of a few x 10~^ Mq yr~^. We also find that the predicted modified 
disk emission is compatible with existing near-infrared (NIR) observations of 
Sgr A* in its quiescent state provided that the disk inclination angle is >87° and 
that the jet extracts more than 75% of the accretion power. 

Subject headings: Galaxy: nucleus - accretion - black hole physics - galaxies: jets 
- MHD 



1. Introduction 



The exceptionally compact radio source Sgr A* is spatially coincident with the dynamical 
centre of the Galaxy, which co ntains a mass M ^ 3.7 x lO^Mpj , deduced from measurements 



of stellar orbital motions (e.g. iGhez et al.ll200d : ISchodel et al.ll2002l : lEisenhauer et al.l 120031 ). 
High- resolution VLBA imaging of Sgr A* at millimetre (mm) wavelengths provides the most 
compelling observational evidence yet for the existence of superma ssive black holes (SMBHs), 
with M now constrained to lie within a radius ~ 1 AU of Sgr A* ( IShen et al.ll2005l ). 



The observed bolometric luminosity of Sgr A* , L ^ 10^^ erg s^^ ^ 2 x lO^^LEdd (where 
Lehh is the Edd ington luminosity), is unusually low, even compared to that of other LLAGN 
(IHo et al.l 119971 ). This can be attributed to a very low mass accretion rate. Ma, or a low 



- 2 - 



radiative efficiency, e = LjM^(? ^ or to a combination of both. The most popular accretion 
models for Sgr A* are Radiatively Inefficient Accretion Flows ( RIAFs), which have attributed 
the low L to a low e and which are base d on the hot ion torus (IShapiro. Lightman fc Eardley 



19761: 



Ichimarul 1977 : Rees et al. 1982 ) and Advection Dominated Accretion Flow (ADAF; 



Narayan fc Yi 



19941 ) models. 



Applications of RIAFs to Sgr A* have evolved con siderably in recent years in response 



to increasing observational constraints (see lYuanll2006l for a review). In the original models 



the inferred mass accretion rate was comparable to the Bondi rate for spherical accretion, 
~ lO~^M0yr~^, and the radiative efficiency was very low, e ^ a few x 10~^, as a result of 
prefere ntial ion viscous heating and negligible electron-ion coupling. The most recent RIAF 



model (lYuan. Quataert &: Narayarul2003l . |2004| ) now takes into consideration a mass outflow 



component as well as energization of the electrons. As a result, the revised RIAF model now 
attributes the low luminosity of Sgr A* to a combination of low Ma and moderately low e. 
The outflow reduces the rate at which hot ionized gas is fed onto the central SMBH to levels 
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requi red to produce the observed levels of radio and sub-mm emission (j Quataert fc Narayan 
19991 ). resulting in a radiative efficiency e ^ 10~^. 



The chanelling of binding energy away from the ions as a result of outflows and electron 
heating (as well as electron-ion coupling) necessarily implies that the ion energy is subvirial. 
There may then be insufficient ion pressure to support a geometric ally-thick structu re and 
consequently, some of the gas can radiatively cool on an inflow time (IRees et al.lll982l ). This 
raises the possibility that a compact, cool accretion disk may form at small radii. To test 
this possibility, the effects of a mass outflow on the continuity, angular momentum and 
energy equations need t o be self-consistently consid ered in the exisiting RIAF model for 



Sgr A * (lYuan et al.ll2003l ). Numerical approaches (e.g. iHawley fc Balbusll2002l : ISharma et al. 
20071 ) have so far been limited by the nonconservative numerical scheme and the neglect of 
radiative cooling. 

Observationally, increasing evidence is emerging to suggest that geometrically-thin disk 
accretion may persi st at or near t he last marginally stable orbit in low- luminosity sources 
jMiller et al.ll2006al jbl: [Maozll2007l : iRvkoff et al.l[2007l : IlIu et"aDl2007h . Yo ung massive stars 
seen orbiting the Galactic Center (GC) and belieyed to have formed in situ (INaYakshin fc Sunvaev 



2006 



2005 : iNavakshin et al.ll2006l : iPaumard et al.ll2006l : lNavakshin]l2006l : iLevin et al. 
2006h suggest that at least some of the hot gas detected by Chandra has condensed into a 



Beloborodov et al 



cool phase. Additionally, new models for stellar wind accretion in Sgr A* indicate that slow 
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winds can r adiatively cool with in a dynamical timescale to produce a cold phase of ac- 
creting gas (ICuadra et al.ll2006l ). However, a strong constraint on the presence of a cold, 
optically-thick disk is the absence of eclipses in the orbit of S2, the best studied star in the 
central S-cluster. This limits the size of a putative accretion disk to < 10^^ cm ~ 3 x lO^rg 
(ICuadra. Navakshin fc Sunvaevll2003h . where = GM/c^ ^ 5 x 10^\M/3.7 x lO^M©) cm 
is the gravitational radius. 

If a compact, cool disk forms near SgrA*, much of its emission would suffer strong 
interstellar extinction. Even so, the mass accretion rate would have to be extremely low 



and the spin axes of the disk and galactic plane would have to be closely a. 
emission from a standard disk to remain below observational limits in the NIR ( 



2003 



Falcke 



2001 



Trippe et al.ll2007h. However, it has been suggested (e.g.lNagar. Wilson &: Falcke 



Falcke. Kording fc Markoffl 12004 : iGallol 120071 : Merloni fc Heinzl 120071 ) that the 



igned for the 
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spectral properties of LLAGN and galactic X-ray binaries in th eir low/hard X-ray stat e 
can be attributed to a jet-dominated mode of accretion (see also iMerloni fc Fabianll2002l ). 
Although jets have not been observed in Sgr A* , various arguments have been made for 



their presence (see e.g. iMarkoff. Bower &: Falckd 120071 . for a comprehensive discussion) 



If jets are present, an acc r etion disk can no l onger be described by the standard model 
( IShakura &: SunyaevI 1 19731 : iNovikov &: Thornd Il973l ) because the radial disk structure is 
modifi ed by the magnetic torque responsible for jet formation (IKuncic &: Bicknelll |2004J . 
2007bl ). Indeed, extraction of accretion power by a magnetized jet resu lts in an accretion 
disk th at is less luminous than a standard disk accreting at the same rate (jKuncic &: Bicknell 
2007al ). That is, the radiative efficiency is lower than that predicted by standard accretion 
disk theory. 

In this paper, we revisit theoretical constraints on the formation of a cool accretion disk 
and we investigate observati onal constraints oii jet-d r iven disk accretion in Sgr A* . U nlike the 
jet-RIAF model for SgrA* jFalcke^Markofj bood : |Yuan. Markoff fc Falckd l2002h and the 
disk-corona model for LLAGN ( Merloni &: Fabian I2OO2I ) . we explicitly model the magnetic 
coupling betwee n the accretion flow a nd the outflow. This model has been successfully 
applied to M87 (jJoUey fc Kuncid 120071 ). The organization of this paper is as follows. In 
Section [2], we show that a RIAF becomes geometrically-thin as a result of diversion of 
binding energy from the inflowing ions. This allows the electrons to radiatively cool and to 
thermalize with the ions on an inflow time, thus forming a cool disk at small radii. We also 
present the relevant equations for relativistic disk accretion modified by MHD stresses and 
we calculate the modified disk fiux radial profile and corresponding modified disk emission 
spectrum using parameters appropriate for SgrA* . In Section[3], we calculate the steady-state 
synchrotron spectrum resulting from a jet magnetically coupled to the underlying accretion 
fiow. We compare the predicted spectra with the observed quiescent spectrum for Sgr A* . 
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A discussion and concluding remarks are given in Sections H] and [5l respectively. 



2. Coupling a Magnetized Jet to a Relativistic Disk 

2.1. Formation of a Cool Disk 

Weakly magnetized, differentially rotating accretion ffows are unstable to MHD turbu- 
lence and electron heating is unavoidable as a result of viscous dissi pation via a turbulent 
cascade as well as resistive dissipation v i a stochastic recorinection (IQuataert &: Gruzinov 
I999I : iBisnovatvi-Kogan fc Lovelacel I2OO0I : ISano fc Inutsukal I2OO1I : Isharma et aliboOTh . As 
the accretion flow becomes increasingly radiatively efficient, the internal energy of the ions 
drops to subvirial levels and the accretion flow geometry deflates from a geometrically-thick 
torus (iRees et al.l Il982l ). If an outflow is also present, then this can further promote the 
collapse to a geometrically-thin structure, with a height-to-radius ratio h/r <^ 1. 

To see this quantitatively, consider the internal energy equation for ions, with specific 
energy Ui, in an axisymmetric, steady and incompressible accretion flow with radial velocity 
Vr in which a fraction 6 of the binding energy extracted by the internal MHD stresses tr^ is 
diverted to the electrons and in which there is mass outflow with velocity Vz'- 
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on 
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Note that the current RIAF model (I Yuan et al.l |2003| ) neglects the outflow term involv- 
ing Vz- Vertically integrating and using the relations Ma = 27rrS(— f^) and dM^/dr = 
—ATTrp{h)vz{h), where is the mass loss rate and where continuity implies dMi^/dr = 
—dM^/dr, gives 

MaMi ) + Ui— = -2TTrTr^ T— (2) 



dr \ / dr or 

where T^,^ is the vertically-averaged stress. This is obtained from conservation of angular 
momentum: 



2'KrTr. 



1 - 



MJr) 



)"•] 


-'-1 







r'Q- 



.dM^ 
dr 



dr 



(3) 



where r^s is the last marginally st able orbit of the inflow. Substituting this into ([2]) to 
ehminate Tr^ and using Ma oc r* (lYuan et al.ll2003l ). where < s < 1 and a keplerian 
angular velocity, Vt = (GM/r^Y^"^, yields the following solution for the ion internal energy: 



Ui{r) 



GM 



(4) 
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where ((r) = 1 — (3/2 + s) ^(r/rms) is a small-r correction factor. This solution can 

be written in terms of the ion temperature and the ion virial temperature T^iri'- 



JL 

T ■ ■ 

-'- vir,i 



r(i-^) 



(5) 



Thus, the ions are subvi riaL The current RIAF model for Sgr A* requires 5>0.5 and nu- 
merical simulations (e.g. iHawley &: Balbud |2002| ) indicate s ~ 0.5 — 1. Taking s = 0.75 for 
example, gives Ti/Tvir,i ^ 0.2. 

This result has important implications for the radiative cooling and electron-ion col- 
lision rates, both of which are very sensitive to the scaleheight ratio h/r ^ T^/T^i^^i of 
an ion-pressure-supported disk. Con sider the inflow timescale, f^^"^ ^ ''"/I'^rl, and the 
bremsstrahlung cooling timescale (e.g. Rybicki fc LightmanI 1979 ). 



^cool ^ 



7 X 10^ V-'Ty2 s 



(6) 



where rig is the electron number density (in units cm ^) and Tg is the electron temperature 
(in K). Using ng ~ M^,/ [Attt ^m.ph\vr\) the ratio of cooling to inflow timescales is 



j-cool 



^inflow 



e h vl 
- 100 



1/2 



(7) 



^Edd 



IS 



|aCsr from the a-disk formalism, 



where m = M^,/ M^dd is the dimensionless mass accretion rate and Msdd 
the Eddington accretion rat^. Using the relation \vr\ 
where Cg ~ {jkTi/ fiiripY^'^ = Vth is the isothermal sound speed, the condition ^cooY^mfiow < ]^ 
required for the electrons to cool before reaching the black hole then implies (see also e.g. 
Rees et al.lEgS^ 



m>2 X 10"^e_2a^2 



h 



1/2 



, . o , (8) 
O.lr / yrrif.c'^ J 

where e_2 = e/10~^ a nd 0^-2 = a/ lO"^. Typical values of a found in numerical simulations 
are a few percent (see lBalbusll2003l ) . Note that m = L/L^dd ~ 2 x 10~^ for Sgr A* , implying 
that hot electrons with initial temperatures Tg < lO^'' K can cool down by the time they reach 
the black hole. 

If radiative cooling by bremsstrahlung emission resulting from electron-ion encounters 
can occur on an inflow time when the ions are subvirial, then it is necessary to also reconsider 



^Note that we define MEdd without assuming a 10% radiative efficiency, as is done in RIAF models, in 
order to keep the dependence on e exphcit throughout the equations. 
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whether the Coulomb colhsion time can be shorter than the inflow time. The electron-ion 
collision time, t*^' ~ (7r/2)-'^/^mp/(mcncCrTclnA)(fcTe/meC^)^''^, can be written in terms of 
accretion parameters as 



7r\i/2 1 e hlvrltTLp r f kT^ 
2 



In A m rrie Vg \ rUpC^ 



3/2 



(9) 



The condition < t '^'^ then requires 



m > 50 e a 



5 X 10"^°e_2a^2 



h 



O.lr 



3/2 



(10) 



where a Coulom b logarithrn InA = 50 has been used. Note that this differs from previous 
calculations (e.g. iRees et al.lll982l ) only in the explicit inclusion of h/r (previously taken as 
unity) and e (previously taken as 0.1). For Sgr A* , with m ^ 2 x 10^^, this implies electrons 
with temperatures Tg < 2 x 10^° K can thermally couple with subvirial ions on an inflow time 
via Coulomb collisions. Note that collisionless wave-particle plasma microinstabilities can 



enhance the electron-ion coupling rate further (IBegelman fc Chiueh 



2000 



MeliaetaL 



200l|). 



1988 



Quataert 



1998 



Gruzinovl 1 19981: Quataert fc Gruzinov 1999 : Blackman 1999 : Bisnovatyi-Kogan fc Lovelace 



As the disk cools, its scaleheight continues to decrease until the gas becomes optically 
thick. Initially, the opacity of the hot gas is dominated by electron scattering. The optical 



depth over the disk height is Tes ~ crTAfa/(47rr/imp 



Ma>7 X 10"^a_2 



1/2 



lOOrg/ VO.lr 



which exceeds unity when 



MeM^yr- 



(11) 



For Sgr A* , this implies a disk can become optically-thick once it collapses to scaleheights 
/;,<afew x 10~^r. By then, however, the disk is sufficiently cool and dense that other 
opacities are more important than electron scattering, so the final disk height need not be 
too small. 



We note that whilst numerical simulations (e.g. ISano &: Inutsukal l200ll : ISharma et al. 



20071 ) demonstrate that preferential heating of electrons resulting from viscous and resistive 
turbulent dissipation occurs in low-Ma accretion flows, they are yet to include radiative 
cooling in the calculations. This is needed in order to verify the above theoretical arguments. 
In the following, we assume that a cool, quasi-thermal disk has formed at small radii as a 
result of electron heating and mass outflows in the hot accretion flow at larger radii. A 
thermally-driven outflow is assumed to be negligible in the cool disk, but we consider the 
effects of a magnetically-driven outflow which we identify as a magnetized jet. 
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2.2. Magnetic Torque on the Disk Surface 



The MHD stresses responsible for allowing accretion to proceed can give rise to a non- 
zero torque that acts o yer the disk surface, in addition to a non-zero torque at the last 
marginally stable orbit (iKuncic fc Bicknellll2004l ). The effect of a non-zero magnetic torque 
acting on the surface of a relativistic disk has not been self-consistently modelled before. 
This torque does work against the disk surface, thereby removing energy from the disk and 
directing it vertically. Because it is associated with the azimuthal-vertical magnetic stress 
tj^ = B^B^ /An at the disk surface (denoted by the '+' superscript), it can be identified with 
the formation of a magnetized corona and/or jet. Thus, the radiative fiux from a relativistic, 
torqued disk can be most generally expressed as (see the Appendix for details) 



Fir) 



3GMM, 



[r"(x) + /.T(^)-/;f(^)] 



(12) 



where x = r /r^ is the dimensionless radius, /^^(x) is the iNovikov fc Thornd (119731 ) rela- 
tivistic correction fac tor, f^f'(x) is a cor rection factor for a nonzero torque (NZT) at the 



inner disk boundary (lAgol &: KrolikI l2000l ) . and f^Tix) is an analogous correction factor for 



a nonzero torque on the disk surface. This is derived in Appendix lA. 21 and is explicitly given 
by 



nzt / 



B 



-ttx dx 



(13) 



where x^s = i^ms/fg is the dimensionless radius of the last marginally stable orbit and the 
dimensionless functions C{x) and B{x) are defined in the Appendix. 

Henceforth, we consider a jet interpretation for the nonzero surf ace torque and assura e 
a simple power law radial profile for the (j)z stress as follows (see also iFreeland et. al.ll2006l ): 



(14) 



The normalization of the nonzero MHD stresses t 



(j>Z 



Xr 



and Tr^{xms) are yet to be de- 
termined. We now achieve this by considering the global constraint imposed by energy 
conservation. 



2.3. The Global Energy Budget 

Global energy conservation requires that the total accretion power Pa niust equal the 
sum of the disk radiative power and the jet power Pj, which we equate to the rate at 
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which work is done against the disk by the nonzero magnetic torque on its surface. Hence, 

Pa = + P-i (15) 

We now write 

Pa = eaMac' = (e"" + Ae)Mac2 

Pj = ejMaC^ (16) 

where ea = e'^'^ + Ae is the total accretion efficiency and ej is the efficiency of energy removal 
to a jet. The efficiencies are defined as follows: 

= I r x-^rix) dx = ir x--^ dx (17) 
^^=\r ^-'W - gg 37rxL7.,(x^s)r, ^^^^ 

^ Jx^s ^ms MaC 

= I --Vr (-) = '"""'f^'"--' /2 (19) 



where 



roc 5/2 

Thus, the nonzero torques can be normalized in terms of the corresponding efficiencies Ae 
and ej and the terms /°|*(a;) and f^f'{x) can be rewritten as follows: 



-9- 






Fig. 1. — Radial flux proflles for the jet-modifled disk model (solid curves) with different 
physical parameters: a is the dimensionless black hole spin parameter; ea is the overall 
accretion efficiency, which includes a fractional contribution of 5% of the Novikov-Thorne 
efficiency resulting from a nonzero torque at the last marginally stable orbit r^^s', and ej/ea 
is the fractional power removed by a surface magnetic torque responsible for jet formation. 
The disk luminosity is — 10~^LEdd ~ 5 x lO^^ergs"^ in all cases. The dotted curves show 
the corresponding radial flux proflles for a standard, non-torqued relativistic disk. 
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2.4. The Modified Disk Flux Radial Profile and Spectrum 



We now substitute the above expressions into the disk flux proflle f|T2|) : 



(25) 



87rr3 [C{x) 3C{x)xy^h SC{x)x^/^ h 



The global energy constraint ( JTSl) yields the following for the disk radiative efficiency: 



where L^dd = 47rGMmpc/crT ~ 1.3 x lO^^Mg erg s~^. Thus, once Ae and ej/ea are specified 
(and note that e^^ is given explicitly by ( fTTl) and depends only on the dimensionless spin 
parameter a), we have determined ed- Then the mass accretion rate is determined by speci- 
fying Ld/L^dd and using fl271) . The fraction Cj/ea of accretion power injected into a jet cannot 
be arbitrarily large, however, since the local disk fiux must remain positive at all radii. So 
for a given set of input parameters, a, Ae/e^"^, and Ld/L^dd, there is a maximum allowable 
value of Cj/ea. 

The radial fiux profiles predicted by our torqued relativistic disk model (c.f. (125|) ) are 
shown in Figure [1] (solid curves). In all cases, the rate of radial decline in the surface 
magnetic stresses is fixed at g = 2.5 (c.f. (IT^ ). The efficiency of energy dissipation by the 
magnetic torque at the inner boundary is set to Ae = 0.05e^'^. The intrinsic luminosity of 
the modified disk is Ld = 10~^LEdd ~ 5 x lO'^^ergs"^ in all cases. The fiux profiles are 
calculated for a range of black hole spin parameters a and fractional jet powers ej/ea, as 
indicated in each plot. Note that as a, and therefore ea, increase, a larger fraction ej/ea of 
the accretion power must be electromagnetically extracted to maintain a constant Ld- The 
mass accretion rates are as follows: (a) Ma ~ 3 x lO~^M0yr~^, (b) Ma ~ 1 x lO~^M0yr~^ 
and (c) Ma ~ 1 x lO^^MQyr^^. The corresponding profiles for a non-torqued disk (dotted 
curves) are calculated using the same values of Ma, but in that case, ej ~ ea because all the 
accretion power is dissipated locally in the disk. 

The effect of the nonzero torque across the disk surface is to do work against the disk, 
thus reducing the disk fiux over the range of radii where the magnetic torque is strongest. 
This counteracts the effect of the nonzero torque at the inner disk boundary, which enchances 
the disk flux near a;ms. The combined effects of these two torques is clearly evident in the 
radial flux proflles in Fig. [T] (solid curves). The resulting radial proflles are substantially 




(26) 



We can also write 




(27) 
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modified from tlieir corresponding zero-torque profiles (Fig. [H dotted lines). It is clear from 
Fig. [1] that the nonzero magnetic torque acting on the disk surface results in a disk radiative 
efficiency ed that is lower than that of a non-torqued disk. 

Figure [2] shows the disk spectra corresponding to the flux profiles in Fig. [H The spectra 
have been calculated assuming local blackbody emission in each annulus of the disk, with 
the local temperature determined by T(r) = [F{r) / aY^^ ^ where F{r) is given by (!25l) and a 
is the Stefan-Boltzmann constant. We have included a cos^d projection factor, where ^d is 
the angle between the disk spin axis and our line of sight. Note that in all cases, the effect of 
the MHD torque on the disk surface is to make the emergent spectrum (solid curves) dimmer 
and redder than that of a non-torqued disk accreting at the same rate (dotted curves). The 
effect is strongest for the high spin case (Fig. ^) because the overall accretion efficiency is 
highest in that case and to maintain the same Ld, the jet must remove a larger fraction of 
the total accretion power. 



3. Jet Synchrotron Emission 



We identify the nonzero magnetic torque across the disk surface with the mechanism 
responsible for extracting accretion power from the disk and injecting it directly into a jet. 
This torque could also produce a magnetized corona l outflow, in which case only a fraction 
of the extracted power may result in a jet (see e.g. iMerloni fc FabianI 120021 ) . Some of the 
magnetic energy is subsequently converted into kinetic energy. We expect a fraction of the 
particles to be accelerated to nonthermal, relativistic energies. Synchrotron radiation by 
relativistic electrons will then con tribute significantly to the radio emission. It has been 
proposed (IFalcke fc MarkoffI |2000| ) that highly efficient acceleration of electrons and syn- 
chrotron emission at the base of the jet may explain the observed sub-mm excess emission 
from Sgr A* . 



3.1. The Jet Model 

We divide the jet into a series of quasi-cylindrical sections of thickness A^, and calcu- 
late the total emission spectrum by summing the contributions from each component. The 
geometry is illustrated in Figure [31 

The jet plasma is assumed to have a magnetic field strength B and to contain nonthermal 
electrons, with Lorentz factors in the range 7min < 7 < 7max5 where 7min and 7max are the 
minimum and maximum electron Lorentz factors respectively. The nonthermal electron 
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Fig. 2. — Predicted disk spectra for the jet-modified disk model. The dimensionless black 
hole spin parameter a is and 9^ is the inclination angle of the disk. All other model parameters 
are the same as those used in the corresponding flux profiles shown in Fig. [H The solid 
curves are the predicted spectra for a relativistic disk torqued at the inner boundary and on 
its surface; the dotted curves are spectra for a non-torqued relativistic disk (i.e. a Novikov- 
Thorne disk). In each plot, the upper two curves are for a disk oriented face-on (6'd = 0) and 
the lower two curves are for a disk oriented nearly egde-on (^^d = 89°). 
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jet axis 




Fig. 3. — Schematic diagram of the jet geometry (not to scale). The jet begins at a height 
Zq and extends to Zj, with half opening angle 0. It is divided into cylindrical sections of 
thickness Az and radius r z(f). The angle between the jet axis and the line of sight is 
and the path length through each cylinder is As. The distance to the galactic centre is 
Dgc ~ 8.5 kpc. 
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energy distribution is given by Nc{'~f) oc 7 ^ where N(, = J Ne{'j)d'j is the total electron 
number density. 

We consider a relativistic jet with bulk Lorentz factor Fj, half-opening angle (p ~ V^j; 
and Doppler factor 

5={T,[l-il-TrY\ose,]}-' , (28) 

where 9j is the angle between our line of sight and the jet axis (see Fig. [3]). We use the 
following simple radiative transfer model to calculate the observed specific luminosity due 
to the net contribution from each jet component (assuming isotropic emission in the source 
rest frame): 



'^obs / ^ '^obs \ / 



(29) 



where AA ^ nrAz sin 6 is the projec ted surface area of each em itting cylinder, S^^^^^ is the 
synchrotron source function (see e.g. iRybicki fc Lightmarull979l for relevant formulas) and 



r 



syn 



I'obs 



6~^kT As 



is the synchrotron optical depth along a path length As ^ 2r/ cosOj, where r zip through 
each section. 

We use an equipartition factor /gq = Ub/Uq to relate the magnetic field energy density 
Ub = B^/8n directly to the relativistic electron energy density Ue = |(7)A^e"^eC^, where (7) 
is the average Lorentz factor. The electron number density and hence the magnetic field 
B decline with jet height z according to 



Ne{z) OC Z' 



B{z) OC z-^ 



The total jet power is 



=F,(1-F-)V^ 



c X 



(Fj - l)iVempc' + -FjiVe(7)mec2 (1 + 2/eq) 



(30) 



where the first term in square brackets refers to the bulk jet kinetic energy and the second 
term refers to the electron kinetic energy and the magnetic energy. This is used to normalize 
the electron number density A^'e,o at the base of the jet once the ratio ej/ea (or equivalently 
Pj/Pa) is specified. 
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3.2. Application to Sgr A* 



Figure H] shows the best-fit broadband spectra predicted by our magnetically coupled 
dis k-jet model plotted against data points for observations of Sgr A* in q uiescence (compiled 
by iNarayan et al.l Il998l : iFalcke et al.l Il998l : iBroderick &: NarayanI |2006| ) . The best-fit disk 
parameters are listed in Table 1. The disk inclination angle is 6^ = 87° in both the zero- 
spin and high-spin cases. We generally find that the jet needs to extract the maximum 
fraction ej/ea of accretion power from the disk to produce an energeti cally significant spectra l 
component in the sub-mm. For the jet spectrum, we used 6j = 55° ( iFalcke fc Markofa 120001 ) 
and we found that a total jet length of just ^ 350rg ^ 6 x 10~^pc was needed to fit 
the lowest radio frequency data points (this corresponds to an angular size of 6' ~ 1.5 mas 
at the distance to Sgr A*). The jet is launched from an initial height zq = O.lrms above 
the disk midplane, where r^s = 6rg for a = and rms = l-5rg for a = 0.99. The electron 
energy power-law index is p = 1.9. The other jet input parameters used are as follows: (a) 



case) Fj = 1.05, /eq = 30, 7, 



30, 7n 



150; and (b) (a = 0.99 case) Fj = 1.9, 



eq 



1 ) Tmin 



50, 7. 



100. The corresponding electron number density and magnetic 



field strength at the base of the jet are: (a) Ae,o ~ 1 x lO^cm"^, Bq ^ 230 G; and (b) 
Ae,o ~ 1 X 10^ cm~^, Bq ^ 430 G. Despite the large difference in Ae,o, the predicted jet 
spectra for the two cases are remarkably similar. This is because the lower Ae^o for the a = 
case is offset by the larger length scale, since we set Zq = O.lrms in both cases and is 
larger by a factor of 4 for the a = case. The location of the synchrotron self-absorption 
turnover is sensitive to the value of Ae^o and we find that zo = O.lrms gives a turnover in the 
sub-mm in both cases, which is suggested by the observations. 

Note that in both the a = and a = 0.99 cases, we find that a relatively narrow range of 
electron Lorentz factors, 7mm $7 $7max5 is required to produce a prominent "bump" feature 
in the sub-mm. The range of electron energies is not sufficiently wide to produce a broadband 
optically-thin synchrotron power-law spectrum, so the spectrum cuts off sharply above the 
maximum critical frequency at the b a se of t he jet, ~ 4 x 10^'y^^^B Hz. A similar requirement 
was also found by IFalcke fc Markoffl (120001 ) for their RIAF-jet model. We find it difficult to 
fit the peak in the sub-mm excess feature at frequencies > 100 GHz. It is likely that one or 
more of our simplifying assumptions (e.g. keeping 7min and 7max constant throughout the 
jet) may need to be relaxed in order to test more thoroughly the jet interpretation of the 
sub-mm excess. However, the observational data are mostly upper limits in this band and 
ther e may be other cont ributions to this spectral component, such as warm dust emission 



le.g. 



Becklin et all I 19821 ). 



The main difference between the high-spin and low-spin disk-jet models is that the 
0.99 model predicts a larger jet kinetic power and a bluer disk spectrum compared to 
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Table 1: Disk parameters used in the coupled disk-jet spectral model shown in Figure HI The 
dimensionless black hole spin is a, Cj/ca is the fraction of accretion power extracted by the 
magnetized jet, ea is the total accretion efficiency, is the disk radiative efficiency. Ma is the 
mass accretion rate in Mq yr~^ and Pj is the jet power in erg s~^. The black hole mass is fixed 
at M = 3.7x 10^ Mq and the total disk luminosity is fixed at = 10~^LEdd ~ 5 x lO^^erg s"\ 



input parameters inferred parameters 

a Cj/ea Ca ed Ma ^_ 

0.00 0.75 0.06 0.02 5 x lO'^ 1 x 10= 
0.99 0.90 0.31 0.03 3 x 10-^ 5 x 10= 




Fig. 4. — Broadband steady-state spectra predicted by the magnetically coupled disk-jet 
model for (a) a zero-spin black hole and (b) a high-spin black hole, compared against dered- 
dened observations of SgrA*in quiescence. The solid line is the total disk-|-jet spectrum 
predicted by our model, the dotted fine is the jet synchrotron spectrum and the dashed 
line is the modified disk spectrum for an i nclination of 6d = 87 °. The diamonds, t rian- 
gles and squares are data po ints compiled by lNarayan et al.l (Il998l ). iFalcke et al.l (Il998l ) and 



Broderick fc NarayanI (120061 ) . respectively. 
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the a = model. The larger jet kinetic power for the a = 0.99 model arises because a 
spinning black hole can extract more accretion power Pa than a non-spinning hole accreting 
at the same rate (and note in Table 1 that the accretion rate for the a = 0.99 model is 
only slightly lower than that for the a = model, whereas the accretion efficiency is more 
than a factor of 4 higher). A bluer disk spectrum is predicted by the a = 0.99 model 
because the accretion disk extends all the way down to a last marginally stable orbit of 
Tins ~ l-5?"g (compared to r^s ^ 6rg for the a = model) and thus, the peak temperature of 
the accretion disk is higher in the high-spin case than in the low-spin case. For an intrinsic 
disk luminosity Ld ~ 5 x lO^^ergs"^, we find ^^d^87° (i.e. the disk is almost edge-on) is 
required for the disk spectrum to fall below the observational upper limit in the near infra-red 
(see Fig. H]). Importantly, this data point places an even stronger observational constraint 
on a standard accretion disk model, which predicts a higher and therefore requires an 
even higher 6^ for the same Ma. This is because a standard disk model predicts that all 
the accretion power is radiated by the disk, whereas our jet-modified disk model takes into 
account magnetic extraction of accretion power by a jet and hence, predicts a dimmer and 
redder disk compared to that of a standard disk accreting at the same rate. 



4. Discussion 

We have demonstrated that the low radiative output from SgrA*is not incompatible 
with a slowly accreting, geometrically thin cool disk. We have suggested that such a disk 
may form at small radii as a result of mass outfiows and electron heating in the hot part 
of the accretion fiow. These effects divert binding energy from the ions, thus reducing the 
scaleheight, allowing radiative cooling and Coulomb collisions to operate on an infiow time. 
Numerical simulations of low- Ma accretion fiows that include radiative cooling are needed 
to test this prediction. We infer cold gas accretion rates of Ma ~ a few x lO~^M0yr~^. 
The disk radiative efficiency, ed ~ 0.02 — 0.03, is lower than the total accretion efficiency ea 
because accretion power is extracted from the disk to form a jet. 

An important observational test of this model is the predicted disk luminosity, which 
i s difficult to constr ain observationally owing to the nearly edge-on orientati on of the disk 



(ITrippe et al.l 120071 ) and strong (~ 30 mag) interstellar extinction (see e.g. ICuadra et al. 
20031 ). We have shown here that an intrinsic disk luminosity Ld ~ 5 x lO^^ergs"^ for 
Sgr A* requires the disk to be inclined at 6'd>87°. Importantly, the constraint on ^d from 
this thin-disk model is not as stringent as that implied by a standard disk because a jet- 
modified disk is dimmer and redder as a result of efficient removal of accretion energy to 
power the jet. The disk luminosity and inclination are constrained by so far only a few 
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flux measurements in the NIR. Clearly, more NIR observations of Sgr A* in quiescence are 
needed to place tighter constraints on the steady-state disk spectrum (see also the discussion 
below on jet X-ray emission). Applications of this model to other more suitable LLAGN 
may thus provide bett er model constraints. So far, we have applied the model to M87 



( jJolley fc Kuncid 120071 ) and have found remarkably good agreement between the predicted 
and observed optical spectra, as well as the predicted jet power and that inferred from jet 
observations on kiloparsec scales. Furthermore, the model suggests that M87 may harbour 
a rapidly spinning black hole accreting at a rate of 10~^MQyr~^. 

We have not attempted here to model the X-ray emission from Sgr A* . However, Comp- 
ton processes could be included in the jet spectral model to calculate the X-ray flux and 
further constrain the physical parameters by comparing with X-ray observations. By consid- 
ering Compton scattering of the disk photons, in particular, we could in principle be able to 
indirectly constrain the intrinsic disk luminosity, L^, and jet-modified spec trum. We note. 



however, that recent simultaneous NIR and X-ray observations of Sgr A* (IHornstein et al. 



20071 ). revealing a constant spectral index during NIR flares and constant X-ray flux, may 
rule out inverse Compton emission by jet electrons off disk photons in favour of synchrotron 
self-comptonisation in the jet as the primary emission mechanism for the X-rays. The con- 
stancy of X-ray emission during NIR flares implies a separate source region for at least some 
of the NIR emission. In our model, emission shortward of 2/im is due primarily to disk emis- 
sion and some of the uncorrelated flaring activity in this band could arise from stochastic 
magnetic reconnection events in the turbulent disk. On the other hand, NIR variability that 



is cor related with variability in the submm and/or X-ray bands (see e.g. lYusef-Zadeh et al. 



20061 ) can arise in our model from prompt acceleration events in the jet that temporarily 



raise the high-energy cutoff in the electron energy distribution. 



5. Conclusions 

We have revisited theoretical constraints on thin-disk accretion in Sgr A* , showing that a 
geometrically-thick, hot, two-temperature accretion flow cannot be sustained in the presence 
of outflows and electron heating. We have revised observational constraints, taking into 
consideration modified disk emission as a result of magnetic coupling to a jet. The magnetic 
torque which drives the jet efficiently extracts accretion power from the disk at small radii. 
This results in an accretion disk that is dimmer than a standard relativistic disk accreting 
at the same rate, so the constraints on disk luminosity and inclination are less stringent. 
For Sgr A* , we infer a mass accretion rate of Ma ~ a few x 1O~^M0 yr"-*^. We find that a 
disk luminosity of Ld ~ 5 x lO^^ergs"^ and an inclination angle >87° are compatible with 
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existing observational constraints. However, more NIR observations of Sgr A* in quiesence 
are needed to provide tighter constraints on the presence of a cool accretion disk. Numerical 
simulations of radiative, low- Ma accretion flows are also warranted. 

The authors wish to thank the referee for comments and suggestions that helped to 
improve the paper considerably. EJDJ acknowledges support from a University of Sydney 
Postgraduate Award. ZK acknowledges support from a University of Sydney Research Grant. 



A. Relativistic Disk Accretion 



The relativistic th eory for steady-state, geo metr ically-thin disk accret ion onto a black 
hole was formulated by lNovikov fc Thornd (119731 ) and lPage fc Thornd ( 1l974j ). In this formal- 
isr n, the radiative disk f l ux can be conveniently expressed as the Newtonian solution derived 
by IShakura fc SunyaevI (Il973l ) multiplied by relativistic correction factors. Specifically, for 
a black hole with mass M and dimensionless spin parameter a, accreting at a rate M^, the 
radiative disk flux can be written as 



F{x) = —M-^M,x-'r{x) 

OTT 



(Al) 



where r = Mx is the Boyer-Lindquist radial coordinate and f^'^{x) is the Novikov-Thorne 
relativistic correction factor, defined by 



r{x) 



A{x) 



where 



X 



1 - 3x"^ + 2ax~^/2 



and 



A{x) 



^ Vms 3a 



y 

ill I 

y 2y Vl/ms 
3(2/2 - 



3(2/1 - a) 



In 



In 



y-yi 

2/2/1(2/1 - 2/2) (2/1 - 2/3) V2/ms - 2/1 
2/ - 2/2 ^ 3(2/3 - a)2 



2/2/2(2/2-2/0(2/2-2/3) V2/ms-2/2/ 2/2/3(2/3-2/0(2/3-2/2) V2/ms-2/3 



with 2/ = x^^"^ , 2/ms = Xyns^^"^ and where 



,1 TT 

2/1 = 2 cos ( - cos « - 2 



A -1 A 
2/2 = 2 cos I - cos « + 3 ) 



In I y 0A21 



1 



2/3 



-2 cos cos ^ a 
(A3) 



are the roots of the equation 2/^ — 32/ + 2a = 0. The last marginally stable orbit in the Kerr 
metric is given by: 



3 + Z2 - sign{a)^{'i - Zi)(3 + Zi + 2Z2 



(A4) 
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where 

Zi = l + (l-a2)V3 [(1 + ^)1/3 + (i_^)i/3] ^ Z2 = sj'ia^^Zl (A5) 



A.l. Correction factor for a nonzero torque at the inner boundary 



Agol fc Krohkl (120001 ) calculated a correction term for the radiative flux from a relativis- 



tic disk torqued at the inner boundary rms- In this case, the solution for the comoving disk 
flux can be generalized to 



F{x) = -M-'M^x-' [r{x) + 
The correction factor for the nonzero torque at x^s is 

Cms 27rXnis ^rt/i (^ms) -^^^ms 



(A6) 



Cix) 



X 



(A7) 



with Cms = C{xnis)j ^ms = ^(^ms), wheie Q = B is the angular frequency of a 

circular orbit at radius r, with B(x) = 1 + ax~^/^ and where 



^rr/)(^ms) d^ 



(A^ 



is the vertical ly integrated radial-azirn uthal magnetic stress at the last marginally stable 
orbit (see also lKuncic &: Bicknellll2004l ). 



A. 2. Correction factor for a nonzero torque on the disk surface 

Here we derive a correction term for the radiative flux of a relativistic disk with a nonzero 
magnetic torque on its surface. The disk is torqued by open fleld lines that do work against 
the surface. We identify this torque as being responsible for the formation of a magnetized 
corona and/or jet and we can write the generalized disk flux as 



(A9) 



where f^f'{x) is the correction factor for th e azimuthal-vertical mag netic stress that gives 
rise to a nonzero torque on the disk surface (IKuncic fc Bicknellll2004j ). 
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We follow the procedure used by lPage fc Thornd (Il974l ). using slightly different notation. 
We introduce the following parameters: 



AirrF 



~m7 



T 



Anrt 



+ z 



(AlO) 



where F is the disk flux measured in the fluid frame and t J ^ is the magnetic stress evaluated 
at the disk surface (IKuncic fc Bicknellll2004l ). 



Conservation of angular momentum can be expressed as 



(All) 



where is the speciflc angular momentum of a circular orbit at radius r. Energy conserva- 
tion gives the relation w = BC^/"^ ^ and substitution into (lAlip then gives 



+ T 



(A12) 



This is a flrst-order differential equation in ^ with a general solution 
C 



^1/2 

^ {Ll - T) dr + 



B 



f-Q 1 



in 



(A13) 



The last term on the right hand side of this equation takes into account nonvanishing mag- 
netic stresses at r^^- Energy conservation and the expression for w in (lAlOp give 



en 



-Bins 2vrrinsT/(rj] 



1/2 
ms 



(A14) 



This term is set to zero in the lPage fc Thornd (119741 ) treatment. 



The term in (]A13p describing the effects of a nonzero suface torque is the one involving 
T. Using {—VL)^r = |-B~^M"^a;~^/^ with ^ and T (see (1A10I) ). the comoving energy flux 
removed from the disk by this magnetic torque is 



+ z 



dr 



(A15) 



from which we deduce that the correction factor f^f'{x) for a nonzero torque on the surface 
of a relativistic disk, as deflned in flAQp . is 



-1 



Anxtt ^ dx 



(A16) 
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